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Abstract: We compute three particle correlations in the Glasma flux tube model of high 
energy heavy ion collisions. We obtain a simple geometrical picture of these correlations; 
when convoluted with final state radial flow, it results in distinct predictions for the near 
side three particle correlation in central heavy ion collisions. 
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1. Introduction 

Recent experiments at RHIC by the STAR |], PHENIX g and PHOBOS § collaborations 
have demonstrated the existence of striking "ridge" like structures in the near side two 
particle correlation spectrum. In near side events with prominent jet like structures, the 
spectrum of associated particles is collimated in the azimuthal separation A<1> relative 
to the jet and shows a nearly constant amplitude in the strength of the pseudo-rapidity 
correlation Ar/ up to Arj ~ 1.5 Q. This collimated correlation persists up to At] ~ 4 |^. 
An important feature of ridge correlations is that the above described structure is also seen 
in two particle correlations without a jet trigger |^. In such events, a sharp rise in the 
amplitude of ridge events is seen ||^] in going from peripheral to central collisions. A number 
of theoretical models have been put forth to explain these ridge correlations 

In Ref. § it was suggested that the ridge can be explained by the radial flow of approx- 
imately boost invariant Glasma flux tubes; gluons are produced with isotropic azimuthal 
distributions in the rest frame of each of these flux tubes. It was shown subsequently 
that STAR results on the centrality dependence of the amplitude of the ridge and of its 

width, for two different energies, could be quantitatively reproduced in this framework 



with only one free parameter |1C]. In this paper, we will compute three particle correlations 
in the Glasma flux tube picture and predict the corresponding near side ridge structure 
of these correlations. Our work is motivated by first studies of three particle near side 



correlations at RHIC [11] 
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The Glasma flux tube picture arises from first principles in the Color Glass Condensate 
(CGC) |l2| , 13| approach to multi-particle production in high energy heavy ion collisions. 
To all orders in perturbation theory, in the leading logarithmic approximation in x, the 
n-gluon correlation spectrum in high energy nucleus-nucleus collisions can be expressed 

as [pira 



(PpL,idyi ■ ■ ■ (Pp^^ndyn 



-{pA, pb) ■ ■ ■ -7) —[pa.Pb)- (1-1) 



d^p±,idyi ' d'^p±^ndyn 

Here pA and ps are the color charge densities of the two nuclei, whose distributions are 
determined by the universal weight functionals W, as described in the CGC framework. 
The weight functionals satisfy the JIMWLK renormalization group equations [16| which 
determine their evolution with the rapidity Y relative to the beam rapidity Ibeam- This 
expression is valid as long as the rapidity gap between two particles is such that \yi — 
Ujl < l/osi where yi and yj are the rapidity of the two observed particles. The rapidity 
Y collectively denotes rapidities in the interval 2/i,---yn- When the rapidity separation 
between any two gluons exceeds this value, additional gluons can be emitted between tagged 
gluons^ . Because our focus here will be on the STAR data, whose pseudo-rapidity coverage 
is limited to Ar] ~ 1.5 units, eq. (LI) is the appropriate expression for our purposes. 

The leading order single particle inclusive distribution, for a fixed distribution of 
sources, is given by 

dNi^o 



d^p^dy. 



iPA.PB) = lim f d^x d'y e'P<^-y^ (9° - iEp){d^y + iEp 



xZ^^iPVxiP) A^{x)ipA,PB)A,{y){pA,PB) . (1.2) 

A 



Here, the gauge fields ^^(x) are solutions of the classical Yang-Mills equations in the 
forward light cone after the nuclear collision. Analytical solutions for these gauge fields are 
known to leading order^ in the sources the full solution, to all orders in pa,Pb, 

requires a numerical computation |21|. In the nuclei, before the collision, the typical range 
of color correlations is the saturation scale Qs, where Qs~^ < Aqcd~^- The saturation 
scale at a given transverse position in the nucleus depends on the two dimensional transverse 
projection of the nuclear matter distribution. In eq. ( |l.lD , it appears in the initial conditions 
for the W functionals; the energy evolution of the saturation scale is determined by the 
renormalization group equations. Because the saturation scales in the two nuclei are the 
only scales in the problem, besides the nuclear radii, the energy and centrality dependence 
of the inclusive observables are a consequence of the properties of the saturation scales in 
the nuclei. 

The expression in eq. (^]^) is remarkable because it suggests that in a single event- 
corresponding to a particular configuration of color sources-the leading contribution is 



^This additional radiation modifies eq. (1.1) significantly For simplicity, we will not quote the full 
expression here but refer the reader to Ref. |l7|]. 

^For interesting recent work in analytical treatments of this problem, see Ref. ||l8|. 
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from n tagged gluons that are produced independently. The coherence in n-particle gluon 
emission is instead generated by averaging over color sources that vary from event to 
event. Because the range of color correlations in the transverse plane is of order 1/Qs, 
this formalism suggests an intuitive picture of coherent multiparticle production as arising 
from fluctuations in particle production from color flux tubes of size 1/Qs from event to 
event. 

The highly coherent matter produced immediately after the nuclear collision is the 



Glasma [22, 23 1. Besides providing the underlying geometrical structure for long range 
rapidity correlations, the Glasma flux tubes also carry topological charge [^4|, which may 
result in observable metastable CP-violating domains |^. While eqs. (|]^) and (^]^) de- 
scribe particle production from the Glasma flux tubes, they do not describe the subsequent 
final state interactions of the produced gluons. If, as widely believed, the produced matter 
thermalizes by final state interactions, these will not significantly alter long range rapidity 
correlations. However, the radial flow of this matter will have a significant effect on the 
observed angular correlations. This is because even particles produced isotropically in a 
given flux tube will be correlated by the radial outward hydrodynamic flow of the flux 



tubes. Elements of these ideas were present in the literature ||2^, 27]-they were proposed in 
the Glasma flux tube picture in Ref. Q and developed further along with analysis of the 
two particle correlations measured by STAR in Ref. |To| ]. 

We will extend the approach outlined above to compute the three particle spectrum 
in nucleus-nucleus collisions in the Glasma flux tube model. The paper is organized as 
follows. In the next section, we briefly review the two particle computation and then 
proceed to discuss the Glasma 3-particle correlation. In section 3, we compute how radial 
flow affects the three particle correlation. In the following section, we discuss our results 
and predictions for experimental three particle correlations. We conclude with a brief 
summary. Technical details of the computation are contained in three appendices. 

2. Computing the Glasma 3-particle correlation 

Before coming to the three particle correlation let us first summarize the result of the two 
particle correction computed in HI . Many of the arguments made here will carry over to the 
three particle case. We will focus in this section on the intrinsic multi-parton correlations 
arising from the physics of high parton densities. The effects of final state interactions, in 
particular the effects of radial flow, will be addressed in the next section. 

For two particles having transverse momentum p_|_ , q_|_ the correlation takes the form 

C (p q) - / '^^^^ \ / dN \ / dN \ 

^ ' ~ \dypd'^p±dyqd'^ci±/ \dypd'^p± / \dygd'^ci^ / ' 

Let us first start by considering a fixed configuration of color sources. Then there will be 
connected and disconnected pieces as shown for example in fig. ||. One might naively think 
the disconnected diagram would simply cancel with the uncorrelated terms in eq. (p.lj). 
This is not the case because the averaging over the color sources brings about non-trivial 
connections between the seemingly disconnected diagrams. 
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Figure 1: Top Figure: A classical diagram which yields a non- vanishing two particle correlation 
after averaging over the color sources. Bottom Figure: A contribution to the correlation function 
associated with a quantum correction to the classical field. In the strong field limit, each source 
goes like l/g while each vertex yields a power of g. One therefore finds that the bottom diagram 
is power suppressed relative to the top diagram by . Analogous arguments hold for the three 
particle case. 



Any diagram that appears to be connected for a fixed distribution of color sources is 
either included in the renormalization group evolution of the sources (if it contains as many 



powers of the rapidity as of Ug) or is else suppressed |14, 15]. The net result is that only the 
classical disconnected graphs contribute; all the effects of quantum evolution are hidden in 
the source distributions. We shall compute these classical contributions and evaluate the 
non-trivial correlations that result when they are averaged over the source distributions. 
In performing the averaging, there will be trivial color correlations which will cancel with 
the subtractions in eqs. (^]^). The final result for the two particle correlation is 

where 

(N^c - 1) m (§^) m (^) 

was estimated analytically in 1^]^. It was found that the particles in a single flux tube, in 
the kinematic range Qs ^ A;^, are produced isotropically in the azimuthal direction around 
the beam axis and two particle correlations are independent of the relative azimuthal angle 
between the tagged gluons. It was then conjectured that this isotropic distribution also 
holds for k±_ < Qs ~ 1 GeV. Recent non-perturbative numerical simulations of Yang-Mills 



''The work as well as a prior version of our work had some erroneous factors of 2 and 27r |28j. We 
have corrected for this in quoting the above value of k,2. 
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equations appear to confirm this conjecture [^] . Tliese numerical simulations also provide 
a more accurate determination of the constant K2. 

The above result has a physically intuitive meaning. If only particles within a flux 
tube are correlated with each other, the ratio of the two particle correlation divided by the 
product of the two single particle inclusive distributions is then simply proportional to the 
ratio of the flux tube area Qg^ to the total system size area S±. If this geometrical picture 
is correct, it should also hold for higher particle correlations. In particular, we will show 
that the three particle correlation can be expressed as 

C3(p,q,l) = [^^^J i^-^^^) ) (2.4) 

where K3 ~ 0.3. The three particle correlation divided by the product of three single 
particle inclusive distributions is simply proportional to square of the ratio of the flux tube 
area over the transverse area of the system. 

We now turn to the explicit calculation of the three particle correlation. We first start 
with the definition of the variance of the three particle multiplicity distribution, for three 
particles having momentum p, q and I. 

c^3(p>q>i) - f.ll ...... ) (2-5) 



dypd'^p±dyqd'^q±dyidH^ 

d^No \ / dN \ I d'^No \ / dN 



dypd'^p±dyqd'^q± / \ dyidHj_ / \ dyqd'^q±dyidH± / \ dypd'^p± , 
d'^N2 \ / dN \ / dN \ / dN \ / dN 



+ 2 



> dypd^p^dyidH^ / \ dyqd'^qj_ / \ dypd'^pj_ / \ dyqd^a^^ / \ dyidHj_ ^ 

In the above expression, the angular brackets (•••), designate an event averaged quantity. 
In our formalism this amounts to averaging over the color sources of the two nuclei. In 
order to compute C3(p,q, 1) we need simply compute the connected contributions to the 
first term in eq. (|2.5| ). 

The classical contribution to C3(p, q, 1) can be computed analytically for Qs ^ P±, 
Just as in the case of the two particle correlations, we anticipate that the result at large 
transverse momentum will demonstrate key features of the correlation that are generic 
and will therefore persist at lower momentum as well, even though the overall numerical 
coefficient in the result will likely differ. This conjecture can be confirmed by numerical 
solutions of the classical Yang-Mills equations; as we discussed previously, ongoing numer- 
ical work for the two particle correlations appear to confirm the conjecture in that case. 
Nevertheless, it is important to note that for large p±, the present formalism may break 
down and other mechanisms such as, for example, beam jet induced rapidity correlations 
may become important. 

The classical amplitude for the production of three gluons having momentum p, q, 1 is 

MTyiUp,^,l) = e^(p)e^'(q)6f (l)/g2/2A-''^(p)yl^''^'(q)yl^-'^"(l) (2.7) 
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where the e's are the polarization vectors of the gluons with the polarization indices A, A' 
and A" and the a, a', a" are color indices for the gauge fields. Following the discussion in 
Ref. 1^ (and references therein), for large transverse momentum the classical gauge fields 
can be expressed as 

J [27ry (p^ - k_L) 

where fabc are the SU(3) structure constants and pi,p2 are the Fourier transforms of the 
color charge densities of the two nuclei. Here is the well-known Lipatov vertex which 
is discussed further in appendix A. Performing the sum over polarizations^, color indices 
and taking the modulus squared of eq. (|2.7| ), dN^ can be expressed as 



dypd^P±dy,Jci^dyidHj) ~ 8(27r)9 ^ ^d^f ifabcfa'defa"fgfahifa'jkfa"lm) J U. ^271^1,^^ 

L,(p,ki^)L"(p,k2x) 



X 



(p^ - ki_L)^ (p± - k2±f 

Lpjci, k3_L)-L^(q, k4_L) 
(q± - k3±)^ (q± - k4_L)^ 

(1± - ksx)' (L - kex)' 

jrbcdefghijklm^^^ q, 1; {k^x }) (2.9) 
where 

X^«(p^ - k2±)pt{q± - k4^)pr(l± - k6±)p^(p± - kix)pi(q± - k3x)/5f (1± - ksx)) • 

(2.10) 

The average in the above expression for corresponds to an average over the color 
sources 

{0) = J [DpiDp2]W[pi]W[p2]0[pi,p2]. (2.11) 
In the MV model for large nuclei [^] , the weight functional is modelled by a local Gaussian 

Wlp].e.,(^-JciW^:^^). (2.12) 

In the above expression p\ is the color charge squared per unit area and can be expressed in 
terms of the saturation scale Qs as Qs ~ 0.6 g'^pA ||3^- For a Gaussian weight functional, 
the source correlator in momentum space is 

(/5-(ki^)/5^(k2^)) = (2^)VA^"''5(ki± -k2^) 

(p'^(ki^)/5^(k2^)) = (2^) Vi 5"" S{k,^ + k2^) (2.13) 



*We work in Feynman gauge: ~ —g^^v 
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Figure 2: Graphical notation used for an emitted gluon of momentum p. 




Figure 3: Example of one of the twenty-seven completely disconnected diagrams. This type of 



contribution to the three particle spectrum cancels in the difference of eq. (2.5) 



At small X, where quantum evolution effects are large, the weight functional is given by 



the solution of the JIMWLK Hamiltonian |16 |. 

We can now finally discuss the evaluation of the three particle correlation. Examining 
the structure of !F in equation ( 2.10I) , we see that there are 15 x 15 = 225 possible quadratic 
combinations of the pi's and p2's. We shall now convince the reader that of these 225 
possible diagrams only 16 diagrams contribute to the intrinsic three particle correlation 
in the limit of Qs P±,Q±,l±- In order to simplify our discussion of the diagrams, 
we introduce, as shown in fig. |2|, a graphical notation for an emitted gluon. Then the 
contraction among the sources can be made by solid lines at the top and bottom of the 
boxes as shown for example in fig. ^ 

As mentioned previously, a large class of diagrams will be disconnected and will cancel 
with the subtractions in eq. (2.5). An example of a completely disconnected diagram is 
shown in fig. |^ There are a total of 27 such diagrams which can be ignored. 

Another class of diagram which can be ignored are those which are self-connected. 
An example is shown in fig. ^. Regardless how one performs the remaining contractions, 
the result will be proportional to d{l± — q±) and are power suppressed. The 5 function 
contribution will be smeared out by final state re-scattering. 

There are still many diagrams remaining. For large , , 1^ the leading contribution 
will come from diagrams with the minimum number of crossings in their contracting lines. 
A similar observation was made for the two particle correlation case. In the three particle 
case there are sixteen diagrams of minimal crossing. Twelve are the diffractive diagram 
as shown in fig. ^ and four are interference diagrams as shown in fig. ^. The diagrams 




Figure 4: Example of a self-connected power suppressed contribution proportional to d{l± — q±) 




Figure 5: One of the twelve diffractive diagrams that contribute to the three particle correlation. 



corresponding to the topologies in fig. ^ and fig. |6| are evaluated in appendices A and 
B respectively. We find that all sixteen diagrams yield the same result; the two sets of 
topologies can be related by simple transformations of the momentum flow in the graphs. 
Multiplying either result in the appendices by 16 we find 



C3(p,q,l) 



S± {g^tiAY^^NliN' 



I) 



327r 



11 



7l /I 



(2.14) 



It will be convenient to express the above result in terms of the inclusive single gluon 
spectrum. This has been shown in p9| , pO| , |3l| , |32| to have the form 

Si. (<7WA^c(iVe'-l). (P^ 



dN 



■In 



(2.15) 



,dypd'^p±/ Att^ p\_ "'\Qs, 

Finally, up to logarithmic corrections, as anticipated at the start of this section, the three 
particle correlation is 

1 / dN \ / dN \ / dN 



C3(p,q,l) = /^s- 



S\Q\ \ dypd'^p± I \ dyqd'^(i± / \ dyidH± 



(2.16) 
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Figure 6: One of the four interference diagrams that contribute to the three particle correlation, 
with 




(2.17) 



At this stage it is difficult to address the theoretical uncertainty in K3 because of the 
logarithmic contributions. The full result requires a numerical computation of the three 
particle correlation from solutions of classical Yang-Mills equations; albeit straightfoward, 
such computations are time consuming. As mentioned previously, because Qs and the 
transverse area are the only scales in the problem, we anticipate that the geometrical 
structure of the result will be robust, with the primary uncertainty being the value of K3. 
Finally, we note that as the number of participants in the collision goes as A'^part ~ S±Q'g, 
our geometrical result predicts that C(p, q, 1) oc l/N^^j.^. 



3. Effect of radial flow 

Before we come to a general discussion on the phenomenological consequences of the above 
calculation, it is important to first consider the effect of transverse flow on the three particle 
correlation. In this section, we will compute the behavior of the three particle correlation 
function when subject to a boost in the transverse radial direction. 

Because the size of the flux tubes in the transverse plane is very small <C S±), 

the particles emitted by a Glasma flux tube experience a common flow velocity. In the local 
rest frame of the flux tube, the three particle correlation at large transverse momentum, 
given by the expression in eq. ( |2.14 ), is independent of both rapidity and azimuthal angle. 



In the absence of flow, we find a flat distribution in Arjpg vs. Arjpi as well as in the azimuthal 
correlation A(j)pq vs. A(j)pi. Even though this result is featureless, it is nonetheless highly 
non-trivial as we shall discuss in section ^. 

The flat Ar] distribution is unaffected by the radial flow of the medium. As we shall now 
describe, additional azimuthal correlations are generated by the flow. As done previously in 
Q, we begin by introducing the particle's rapidity in the direction of transverse flow, Cp,q,i = 
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Figure 7: Schematic picture of three correlated particles emitted from a flux tube of transverse 
size and radial velocity Vr- The transverse size of the collision region is defined as S±. 



— In (tan((/)p_g^i/2)), where the particle's azimuthal angle, </>p,g,i is defined with respect to 
the radial vector pointing towards the emission point-the location of the flux tube in the 
transverse plane. Since all the particles are localized within 1/Qs of each other, we expect 
the same radial boost for all three particles. Fig. |^ shows the setup. 

Because the three particle correlation as computed in the flux tube rest frame is in- 
dependent of the azimuthal angle, the effect of the radial boost enters at the level of the 
Jacobian obtained when one expresses the angular distribution in terms of rapidity vari- 
ables 

^3 ~ ^aII^am = ^^^^ ^^^^ ^^^^ ^fJ^liA ■ ^) 
d(t)pd(t)qd(t)i dCpdCqdQi 

Boosting this expression by an amount Cb = tanh~^ Vr, where Vr is the transverse velocity, 
we find 

^ cosh Cp cosh Cg cosh Ci ^2 2) 

cosh(Cp - Cb) cosh(Cg - Cb) cosh(0 - Cb) ' 

where we use a tilde to denote the boosted quantity. The quantity C3 is evaluated in 
the local rest frame and was given in eq. ( |2 . 1 4| ) -expressed in terms of the single particle 
distribution, it is given by eq. ( p.l6 ). There have been simulations performed by a number 



of experimental groups |3^, A first measurement of three particle azimuthal 

correlations has recently been reported by the STAR collaboration |^^. We expect the 
relevant quantity to be the three particle cumulant, C3 divided by the cubic power of the 
single particle distribution. In addition, it is much more relevant to discuss the correlation 
between the relative angles of the particles. 

We now write the three particle correlation normalized by the cubic power of single 
particle distribution expressed in terms of the relative angles /S.(j)pq = (j)p — (j)q and iS.(f)pi = 
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Vr=0.5 Vr-0.8 




Figure 8: (Color online) Contour plot of A{A(j)pq, A(f>pi) for Vr — 0.5 (left) and Vr — 0.8 (right). 
This plot only depicts the azimuthal structure of ridge like three particle near side correlations and 
does not take into account near side- away side correlations arising from jet fragmentation. 



4>p - (pi- 

The explicit expression for A{A(f)pq, A(/)p;, Cs) is given in eq. ( p. 4]) of Appendix C. It is an 
analytical function of the relative angles and the transverse velocity. It is normalized such 
that, for (b = 0, we have A{A4)pq, Acppi) = 1. 

The function A is plotted in figure |8| for radial flow values Vr = 0.5 and 0.8. The 
radial flow collimates the angular distribution of particles in the direction of the flow. 
This is clearly seen in our results; we find that maximum amplitude always occurs at 
'^(ppq = ^(ppi = 0-all three particles emitted colinearly. This feature of our result is 
consistent with first experimental observations of three particle azimuthal correlations pTj . 
As a corollary to this result, regardless of the magnitude of the boost, we find a minimum 
for A(f)pq = 2tt/3 and A(j)pi = 47r/3 corresponding to the three particles being emitted 
with the furthest possible angular separation. As the boost velocity is increased, the ratio 
between the maximum and minimum correlation increases. To observe these, the strength 
of the away side correlation coming from other "jet-like" mechanisms have to be subtracted. 



4. Discussion 

In the previous sections, we described the consequences of a novel mechanism for multi- 
particle production for long range three particle correlations in heavy ion collisions. The 
underlying picture is quite simple. Boost invariant Glasma flux tubes of transverse size 
1/Qs are formed in heavy ion collisions. They give the leading contribution to particle 
production. (For large rapidity separations, there can be significant violations of boost 
invariance-these too can be computed in the Glasma formalism [jl^.) Particles are pro- 
duced with equal probability along the length of the flux tube and decay isotropically. The 
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Figure 9: (Color online) Preliminary results from STAR of the background subtracted three 
particle pseudo-rapidity correlation in (left to right) d+Au, 40-80% Au-f Au and 0-12% Au-fAu 
collisions. The p± cuts for the trigger and associated particles are 3 < pj_ < 10 GeV and 1 < < 3 
GeV respectively. 



flux tubes flow radially outwards due to the strong pressure gradients generated in the 
collision. While this does not significantly affect long range rapidity correlations, the az- 
imuthal distributions are significantly altered because the emitted particles are collimated 
in the direction of the radial flow. In Refs. IC], this picture was shown to provide a 
quantitative description of STAR data on the near side ridge. Three particle correlations 
provide an additional measure to test this picture. They predict that the strength of the 
correlation is proportional to l/N^^^^. 

In the above classical computation, the resulting three particle correlation has equal 
strength at all rapidity. The above calculation therefore predicts a structure-less correlation 
function of finite amplitude in Arjpg vs. Arjpi. Preliminary observations of such an effect 
have already been made |l^] and shown in fig. A near-side ridge phenomenon is seen 
in the three particle correlation as one goes to more central collisions. In the right most 
plot of fig. ^ a clear jet like structure is seen sitting on top of a structureless plateau that 
extends as far as 1.5 units in rapidity. 

It is the transverse boost that yields non-trivial angular correlations. The PHENIX 



collaboration [^8[ and the STAR collaboration |3£] have extracted an average transverse 
velocity as a function of A'part from blast wave fits to data. We use the PHENIX result as 
input into our boosted flux tube model. We found in the previous section that a maximum 
correlation occurs for Acppq = A(j)pi = while a minimum occurs at A(/)pq = 2tt/2> and 
Acfipi = 47r/3. It is therefore instructive to plot the ratio 6*3(0, 0)/C3(27r/3, 47r/3) versus 
-^part as shown in fig. |l^. We have chosen this ratio because it eliminates the explicit 
dependence on Ug, K3 and Np^^^^. Therefore, the only remaining dependence is on the 
radial flow C,b as a function of centrality. We expect to see this behavior of the amplitude 
irrespective of the rapidity gap between particles. It can therefore be used as a test of our 
model. Alternatively, by fitting the ratio of the three particle correlation, one can extract 
an independent determination of the transverse flow of the system. We should also note 
that the recent experimental data on three particle azimuthal correlations |37] does not 
show "cross-talk" term that would show up as horizontal and vertical strips in the contour 
plot of fig. (^). This result is also consistent with the Glasma flux tube picture because, in 
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Figure 10: Ratio of the maximum correlation C-i{/S.(j)pq = 0, A(/>p/ = 0) versus the minimum located 
at C^{A4)pq = 27r/3, A(j)pi = 47r/3) versus centrality shown as the dashed blue curve. The solid red 
curve shows the ratio at A(j)pq = A(j)pi = tt for comparison. 

the latter, azimuthal correlations are produced by angular collimation of the ridge particles 
by radial flow-this does not permit "cross-talk" terms. Such terms are instead expected 
in scenarios where the angular correlation is sensitive to a ridge formed by jet-medium 
interactions. 



5. Conclusion 

We computed three particle correlations emanating from Glasma flux tubes formed in the 
collision of two sheets of Colored Glass Condensate. A simple geometric interpretation of 
the strength of the correlations is found; the correlation is a number of order unity times 
the square of the ratio of the flux tube size over the system size, or equivalently, it is 
inversely proportional to the square of the number of participants. The correlation is flat 
in Arjpq, Arjpi. Non-trivial azimuthal correlations develop due to the radial transverse flow 
of particles produced by the flux tubes. Finally, we predicted the relative strength of the 
correlation as a function of centrality using a blast wave model. This ratio removes the 
uncertainty from the overall non-perturbative coefficient in our result. It therefore lends 
itself as a method to independently extract the average radial flow of the system. Further 
improvements to our model computation include non-perturbative solutions of the classical 
Yang-Mills equations to compute the three particle correlations, a better hydrodynamic 
treatment of radial flow and hadronization effects, and finally, proper treatment of the 
quantum effects that modify the leading order boost invariant treatment of multiparticle 
production in heavy ion collisions. In particular, combining the Glasma Hux tube picture 
with hydrodynamical simulations is very important for understanding the transverse mo- 
mentum dependence of two particle near side correlations. An interesting first attempt 
combining the Glasma flux tube picture with hydrodynamical evolution has appeared re- 



cently 1^. 
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A. Computation of diffractive contribution 

In this appendix, we shall explicitly compute the diagram in fig. |5|. After performing the 
contraction among the sources of this graph we are left with 

jr(l) = {2'K)'^'^f?^6ci6ekSgm^fdSbjShlS{k;s± - kQJ_)6{k3± - k4^)5(ki_L - k2_L) 

X <5(1_L - k5_L + q± - k3^)5(p_L - ksx + 1± - ke±)5{p± - ki± -<i^ + ki±) (A.l) 



Upon substituting the above expression into eqn. ( ^^ ) we find^ 

(A.2) 

We should note that in order to get the above expression into a form symmetric with 
respect to p^ , q_L and 1^ we had to do a trivial shift in integration variables (ki_L — > 
ki^ -|- p_|_). We now make use of the scalar product of two Lipatov vectors 

L,(p, kx)L"(p, L) = [5'^'^ + e'^e'^] k^ {p± - k^Y ^ (p^ - l^r (A.3) 

and after noting that any angular dependence present in numerator will cancel we are left 
with 

-m^^'^ A^c(A^e-lW^±y (k,^-pJ4(k,^-q^)4(k,^-L)4 

(A.4) 

As we are interested in the limit when p_l, k_L, 1_l S> Qs we keep the term with the fewest 
powers of kij_ in the denominator. In this limit the above expression becomes 



^In evaluating the expression below we have used the convienent property that L^{p±,p± — kix) = 
i'(px,kix) 
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In order to evaluate the above integral we must introduce an infrared cut-off k-^^-^ ~ Qs- 
Qs acts as a natural cut-off as it signifies the onset of non-linear contributions to the 
infrared gluon distributions in the CGC. The final result is 

We should note that sub-leading contributions in individual graphs may have significant 
angular correlations. Our conjecture about the generality of the geometrical picture that 
emerges at high transverse momenta relies on these angular dependences cancelling among 
the (many) sub- leading graphs, as well as being smeared out by non- linear rescattering 
corrections. 

B. Computation of the interference contribution 

We shall here explicitly compute the diagram of fig. ^. After performing the contraction 
among the sources of this graph we are left with 

X 6{l± - k5_L - p± + k2±)6{q± - k^^ k6_L)5(p_L - ki_L - q± + k4_L) . (B.l) 

Substituting the above expression into eqn. (p.Op, we find 



8(2^)11 V- ; -cv-c ^^^A-^j j^i2^ (kix-p±)4(kix-qx)4(kix-U)^ 

(B.2) 

which is the same as expression ( [A.2| ). The subsequent steps are identical. 
C. Angular Integrations 

In this appendix, we work out the angular integrations necessary for evaluating the three 
particle correlation after a boost in the transverse direction. 

First, let us define the three particle distribution and the cubic power of the single 
particle distribution in terms of relative angles. 



C3(A0pg, A(ppl) = / C3{(pp, (j)q, (j)l)6{A(j)pg - (l)p + (j)q)5{A(j)pl - (t)p + (pl)d<j)pd(j)qd(j)l 



ClClCl{A4>pq, A(f>pi) = J Cl{(f)p)Cl{4>q)Cl{4>l)6{A(t)pq -(pp + 4>q)6{A(j)pi -(t)p + 4)l)d(t)pd(j)qd(j)l 

(C.l) 

Now, let us evaluate C3(A(/)pg, A(/)p;). The azimuthal dependence is dictated by the 
Jacobian in eq. (^I^). Upon substitution into the previous expression for C3{A(j)pq, A(j)pi), 
we are left with the following integrals 

cosh C,p cosh C,q cosh C,i 



C3{A(j)pq,A(j)pl) =C3 




cosh(Cp - Cb) cosh{Cq - Cb) cosh(0 - Cb) 

x6{A(l)pq - 4>p + (pg)S{A(j)pl - (l)p + 4>l)d4>pd(f)qd(f)l 

(C.2) 
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where C3 is the local rest frame quantity computed in eqn. or ^J^. The above 

integrals can be done analytically and we can express the result as 



CaCA^pg, AcPpi) = 27rC3^(A(^pg, AcPpiXs) (C.3) 

where 

A{A(l)pq,A(ppi,CB) 

= [25 + 36 cosh(2CB) + 3 cosh(4(^B) — 8(cos A(f)pq + cos{A(ppq — A(j)pi) + cos A(f)pi) sinh^ C,b\ 
X [3 + cosh(2CB) - 2 cos A(j)pq sinh^ Cb]~^ 
X [3 + cosh(2CB) - 2 cos Acjypi sinh^ QbV^ 

X [3 + cosh(2CB) - 2 cos(A(/>pg - Ac/^pj) sinh^ Cb]"^ (C.4) 

The angular integrals in the expression for CiCiCi{A(f)pq, A(f)pi) result in an overall nor- 
malization only, 

C.C.C,..„,A.,,...(^)(^)(^) ,C,5, 



The final result is 

Cz{A4)pq, A4)pl) _ K3 



■A{A(^pq,Ac^pi,CB) (C.6) 



CiCiCiiA<l)pq,Acl)pi) SlQj. 
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